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We report the discovery of a unique object, BD+48 740, a lithium overabundant 
giant with A(Li)=2.33 + 0.04 (where A(Li)=log n Li /nH+12), that exhibits radial veloc- 
ity (RV) variations consistent with a 1.6 Mj companion in a highly eccentric, e=0.67 + 
0.17 and extended, a=1.89 AU (P=771 d), orbit. The high eccentricity of the planet is 
(Vj uncommon among planetary systems orbiting evolved stars and so is the high lithium 

abundance in a giant star. The ingestion by the star of a putative second planet in 
the system originally in a closer orbit, could possibly allow for a single explanation 
to these two exceptional facts. If the planet candidate is confirmed by future RV ob- 
servations, it might represent the first example of the remnant of a multiple planetary 
system possibly affected by stellar evolution. 
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Subject headings: Stars: fundamental parameters — Stars: atmospheres — Stars: 
late-type — Stars: individual (BD+48 740) — Planet-star interactions — Planets and 
CO ; satellites: detection 
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1. Introduction 



There has been a growing number of exoplanets detected around post-main sequence (MS) 
^ | stars. At present, about 50 red giants (RGs) are known to host planetary or brown dwarf-mass 
companions. Recent discoveries also show evidence of multi-planet system s (cf. HD 102272 b, c) 



or multi - brown dwarf systems (BD +20 2457 b,c -iNiedzielski et al.ll2009f) . Although hot jupiter 



type planets exist around sub-giants ([Johnson et aU 2010). the more evolved giants show a paucity 
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of short period and eccentric planets (|Johnson et al.ll2007|) . This is most likely the r esult of tidal 
disruption and/or planet engulfment during the RG phase (e.g lVillaver & Liv io 2009). 



Even rarer is the number of lithium-rich RG stars. According to the standard evolution the- 
ory, when a solar-type star leaves the MS and becomes a RG, its lithium abundance should drop 
from A(Li) ~ 3.3 to about a 1.5-level. In fact, the observed upper giant branch limit is A(Li) < 



0.5 (ILind et alJl2009bh and only a few percent of giants have been observed to have A(Li) > 1.5 
(IKumar et alj (120 1 lh and references therein), with some exhibiting A(Li) ~ 3.3, the value expected 
for a protostellar disk rather than an evolved star. 

In this paper, we report the discovery of a planet in a highly eccentric, long-period orbit around 
a lithium-overabundant giant star. 



2. Observations and data analysis 



BD+48 740 is one of the targets of the Penn State - Toruh Planet Search (PTPS), which is 
devoted to the detection and characterization of planetary systems around stars more evolved than 
the Sun. 

High resoluti on optical spectra dis cussed in this paper were collected with the Hobby-Eberly 
Telescope (HET) (IRamsey et alj|l998l) . The telescope was equipped with the High Resolution 
Spectrograph (HRS) ( Tull 19981) which was fed with a 2 arcsec fiber, working at the R = 60000 
resolution. For the basic data reduction, standard IRAfQ tasks and scripts were used. The signal 

to noise ratio was typically better than 200 250 per resolution element at 5900 A. For the 

Li abundance analysis only one order of the red spectra containing the 7 Li 6708A line was used. 
We have carefully checked all the HET/HRS flat-field spectra for an occasional contamination 
caused by a nearby feature that, depending on the actual RV of the star, may mimic the Li line and 
influence the abundance calculation. All such flat-field spectra were excluded from the analysis. 



3. Results 



BD+48 740 (TYC 3304-90-1) is a V=8.69, B-V = 1.252±0.032, n = 1.36+1.13, K2 giant 
in Perseus. Its basis atmospheric parameters, as well as mass, radius and age were determined by 
Zielihski et al.l (|201ll) and are presented in Table [TJ The stellar rotation velocity was estimated by 
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means of the lFekell (I1997|) method and from the SME spectrum modeling. The estimated range of 
values is given in Tabled] 



3.1. Radial velocities 



BD+48 740 was observed by the PTPS survey between Jan 12, 2005 and March 5, 2012. 
Precise relative RVs were measur ed at 15 epochs over 2229 days using the standard I 2 cell calibra- 
tion technique (|Butler et al.lll996l) . Doppler shifts were derived from least-squares fits of template 
spectra to stellar spectra with the imprinted I 2 absorption lines. 

The measured RV s revealed variations over ~ 90 m s _1 . With the estimated amplitude of 
solar-type oscillations (|Kjeldsen & Bedding! 1 1995|) of ~ 8 ms" 1 and the average RV uncertainty 



of 5.6 m s~ ! , the observed RV amplitude was ~ 10 times larger than the expected observational 
uncertainties. Consequently, we found it justified to assume that the observed RV signal could be 
due to an orbiting planet. 

Several absolute RV determinations for BD+48 740 exist in the literature. Three out of five 



measurements presented in lAbtl (|1970|) agree with our determination, while the other two sug- 
gest the amplitude of 61 kms -1 . The reason for such a scatter in the RV values is probably the 
very strong IS N a doublet at AA 5890, 5895 A, which probably originates in the local IS bubble 
([Welsh et al.ll2010() . and is separated from the stellar lines by 46 km s" 1 . When observed at a res- 

o 

olution of 37 A/mm, the blend shows two peaks wi th a minimum, whose position depends on the 
actual atmospheric conditions. iFamaey et al.l (|2005|) found BD+48 740 to show no RV variations, 
but they computed a different RV for the star. The reason for this discrepancy is not clear. 



3.2. Radial velocities modeling 

The initial, semi-global search of the orbital parameter space to model the observed RV vari- 
ations was carried out with the PIKAIA code (ICharbonneaul 1 1995|) . This was followed by the 
Levenberg-Marquardt, least-squares modeling of the RVs in terms of the standard, six-parameter 
Keplerian orbit within a narrow range selected by the genetic algorithm. All the available measure- 
ments were used and no stellar jitter was added to data. The resulting orbital solution is presented 
in Tableland in FigureQ] Uncertainties in the best-fit orbital parameters (Table© were estimated 
from the parameter covariance matrix. 

To check the reliability of the preliminary orbit, 1000 trials of scrambled radial velocities were 
performed, which resulted in 1 12 solutions with^ 2 lower than the one presented here. Therefore, 
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Table 1: Basic atmospheric and stellar parameters of BD+48 740. 



Parameter 




4534 +24 




2.48 +0.12 


Fe/H 


-0.13 + 0.06 


v m;cro [km s" 1 ] 


1.58 


V macro [km S" 1 ] 


3.2 + 0.5 


v roi sin / [km s _1 ] 


3-3.5 


Vrad [km s _1 ] 


-46.32 


Log L/L 


1.69 +0.20 


M/M 


1.5 + 0.3 


R/R 


11.4 + 0.7 


Prot [d] 


165-192 



Table 2: Orbital solution for BD+48 740 b. 



Parameters 



P [days] 


771.3+ 7.4 


To 


55925.32+6.8 


e 


0.67 + 0.17 


K [m s 1 ] 


36.3 + 9.1 


CO 


140.5+ 47.3 


m 2 sm(i)[M Jup ] 


1.6 


a [Am 


1.89 


x 2 


4.86 


ctrv [m s l ] 


9.26 



-5- 




the corresponding false alert probability for our best-fit model is FAP=1 1 .2%. 

The resulting post-fit residuals are consistent with stellar jitter interpreted as under-sampled 
p-mode oscillations. No correlation between the observed radial velocity and the line bisector 
variations exists (r=0.166) and the correlation between radial velocities and line profile curvature 
(r=0.488) is statistically insignificant. 

The fact that the orbital period is much longer than the estimated rotation period, and the high 
orbital eccentricity make it unlikely that the RV variations are generated by a spot rotating with the 
star. 

We have also searched 3029 WASP photometric measurements of BD+48 740 collected over 
1475 days for any periodic variations and found none within uncertainties. Furthermore, 134 
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epochs of Hipparcos photometry of the star reveal its constant brightness of H p =8.8543 + 0.0017 
mag. 

Therefore, we conclude that the observed RV variations reflect a Keplerian motion of the star 
around the system's barycenter. 



3.3. Li abundance 



The observed spectra were analyzed with the Spectroscopy Made Easy package (SME. IValenti & Piskunov 



199a). SME assumes LTE and parallel plane geometry, but it ignores magnetic fields, molecular 
lines and mass loss. The synthetic spectrum synthesis was based on Kurucz's models of stellar 
atmospheres. 



l,7il(l 




6705 6710 

Wavelength [A] 



6715 



Fig. 2.— Observed HET/HRS spectra (continuous line) of BD+48 740 and BD+42 2315, a low 
lithium giant with A(Li)=0.23, in a narrow range centered at Li I 6708 A. The final SME fit to the 
spectrum of BD+48 740 for the Li line is presented in the inset with the dashed line. 



As an input, SME requires a set of lines from the Vienna Atomic Line Database (VALD, 
Kupka et alJI 19990 identified in the spectrum, the stellar data (including RV), and the instrumental 
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profile. The observed spectrum was considered as a model constraint. Although SME requires 
the M/H ratio, the Fe/H was used as an approximation for the metallicity. The 6695-6725 A 
range was modeled, which covers the 7 Li line at 6708 A as well as several lines of Al, Ti, Si and 
Ca - Figf2J SME was allowed to fit RV, rotation velocity and macroturbulence velocity as well. 
To increase the precision, SME fitting was performed on all 16 spectra obtained with the HET, 
yielding an estimate of the abundance uncertainty (rms). A comprehensive analysis of the lithium 
abundances also requires a correction for non-LTE effects, especially in the c ase of the lithium line 
at 6708 A considered here, which is affected by several non-LTE processes (|Carlsson et al] (119941) 
and references th erein). With SME LTE A(Li)=2.24, and the non-LTE correction of 0.087 from 



Lind et alj (I2009a|) . the final non-LTE value of A(Li)=2.33 + 0.04 (rms) was obtained. This value 



places BD+48 740 on the short list of lithium overabundant giants. 



4. Discussion 



The effective temperature (known quite precisely) allows the placement of BD+48 740 close 
to the stable helium burning part of the evolutionary track of a 1.5 M Q , solar metallicity star of 
Girardi et al.l(|2000|) . Although a relatively small spread in the evolutionary tracks in this region of 
the HR diagram, together with the uncertainty in the luminosity, does not allow an unambiguous 
determination of BD+48 740 evolutionary status, the location of the star on the HR diagram favors 
the interpretation in which the star is evolving toward the Red Giant Branch (RBG) tip. However, 
we cannot completely exclude scenarios in which the star is descending from the giant branch 
toward the horizontal branch (HB) or ascending toward the AGB. Unfortunately our HET/HRS 
spectra do not allow for the determination of the the 12 C/ 13 C carbon isotopic ratio to constrain the 
evolutionary status better. 

Any interpretation of the current state of BD+48 740 past the RGB tip would imply that 
the star had already gone through the helium flash and extended envelope stage. As a result, the 
observed pla net at 1.9 AU would h ave been disrupted by tidal capture or at least have its orbit 
circularized dVillaver & Livioll2009|) . Given the star's evolutionary status and the orbital distance 
of the planet, orbital circularization is not expected for another 0.5 Gyr. 



4.1. Eccentricity 



The observed orbital eccent ricity may be overestimated due to the small number of data points 
available (IZakamska et al.ll201 1|) or it could result from not resolvi ng two components in 2:1 cir- 
cular or near-circular resonant orbits (lAnglada-Escude et al.ll2010|) . If taken as a true value, the 
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observed eccentricity of BD+48 740 b at 1.89 AU, although high, fits into the observed distribu- 
tion of the RV planets as a function of semi-major axis (e.g. iKane et al.ll2012|) . This distribution 
begins to significantly diverge from p urely circular or bits beyond 0.04 AU, and by 0. 1 AU it reaches 
the eccentricity range of 0.0 - 0.5 (|Kane et al]|2012[) . It is interesting to note that the distribution 
of eccentricities in multiple planet systems differs from that of single planets (|Wright et al.l 120091) 
in that it displays lower values. Moreover, when plotted as a function of log (g) of the host star, the 
distribution of orbital eccentricities of planets orbiting giants narrows, showing gen erally mildly 



eccen tric orbits, typically with e<0.45 with only one exception of HD 1690 (e=0.64) (|Moutou et al. 
201 lb . 



Since tidal interaction of the planet with BD+48 740 is still not expected at this point and 
the observed eccentricity is not unusual among planets orbiting MS stars, we could argue that the 
planet eccentricity is primordial, that is, it has not been modified as the star left the MS. However, 
if we consider the fact that the eccentricity is very high for a giant star companion we then need 
to find a mechanism operating only during the post-MS evolution of the star that resulted in such 
eccentricity. 

In principle, post-MS evolution can trigger instabili ty in a multiple system tha t was stable 
during the stellar MS evolution due to stellar mass-loss (|Debes & Sigurdssonll2002|) . However, 
RGB stars are not expected to be losing mass in large quantities and therefore the stability bound- 
ary should not move considerably during this phase. Moreover, a non-adiabatic approximation 
is neither expected to o perate at the current planet orbital distance and stellar mass-loss regime 
(|Villaver & Livioll2007|) nor to excite the eccentricity of the planet. 



It is generally accepted t hat large eccentricities, at any evolutionary stage, have to be excited 



Ford & Rasio 



2008) . Of t he several mechanisms that have been proposed, we find the 



(see e.g. 

Kozai effect dTakeda & Rasioll2005l) very unlikely to take place in BD+48 740 b, since no binary 
companion has been detected with a good degree of confidence. The nearest object to BD+48 
740 is 2MASS 02425927+4855549, a J=13.288 mag star which is separated by 12.3 arcsec, corre- 
spond ing to ~6840 AU (assuming the distance to BD+48 740 d=556 + 129 pc from Famaey et al 



(|2005|) ). The JHK photometry available for the star suggests however, that the 2MASS object is a 
hot, distant field star. Furthermore, the star is not a member of the nearby h andx Per open cluster, 
and its proper motion show that it is not a high velocity star, thus making it very unlikely th at a 
stellar encounter would be responsible for exciting the eccentricity (|Laughlin & Adamsl ll998). 



We are left then with a planet-planet sc attering event as a p lausible mechanism to excite the 
eccentricity of this evolved system (see e.g. iFord & Rasioll2008l) . either during the MS or during 
the post-MS evolution of the star. A viable scenario could be one in which the the BD+48 740 
planetary system contained originally more than one planet. In turn of a planet-planet scattering 
event, the innermost (unobserved) planet was brought into the stellar surface due to either a direct 
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hit or tidal interaction with t he evolving star, while the orbit of the outer plane t was excited to the 
observed, high eccentricity (|Marzari & Weidenschillin g 2002 ; lMardlingll20 1 0|) . 



The derived stellar parameters of BD+48 740 imply a stellar radius of R/Rq = 11.4 + 0.7 
(or 0.05 AU). A small number of con firmed exoplanets around MS stars orbit within this orbital 
distance of their host stars (see e.g. ICollier Cameron et al.l |2007|) . It is very likely that these 
planets, at such small orbital separations, will tidally interact with their host already while on the 



MS w ith the end result of the spiral-in of the planet into the stellar envelope (see e.g. LTackson et al 



2009). More substantially, for Sun like stars, the number of Jovian- type planets in the hot jupiter 
category, orbiting within 0.1 AU, is about 1% ([Mayor et al.l 120 1 It is expected, and indeed 
confirmed by the occurrence f raction of hot jupiters around sub-giant stars with masses M* > 
1.45 M Q (| Johnson et al.1 12010|) . that a large majority of these systems, with the exception of a 
few extreme cases, will survive the MS evolution. However, after the star leaves the MS, and 
especially as it ascends the RGB, the tidal influence becomes important and extends far beyond the 
stellar radius for these almost fully convective stars (|Villaver & Livioll2009l) . 



The BD+48 740 planeta ry system might h ave originally resembled HAT P-13 b, c (IBakos et al 



20091) or HD 217107 b, c (Fischer et all 1 19991) . These are the systems around MS stars, both of 
which have a close-in planet in a nearly circular orbit (a hot jupiter) and another, more distant 
companion in an eccentric orbit. A possible dynam ical evolution le ading to such configuration, 
involving planet-planet scattering, was presented by iMardlingl (|2010f) . The observed enhanced Li 
abundance suggests, however, that in the case of BD+48 740 the planet-planet scattering event 
operated quite recently, during the post-MS evolution of the star. 



4.2. Lithium overabundance 

It is difficult to identify a mechanism that might have preserved the primordial Li-abundance 
of BD+48 740 given that, as soon as the the star moved past the sub-giant branc h, convection 



shoul d have reduced its Li abundance by at least an order of magnitude (see e.g. iGratton et al. 



2000). 



Thus, either BD+48 740 has experienced a Li-production in its interior, or, it has suffered 
external pollution from a companion. We can rule out the presence of a stellar companion that 
might have transfered Li-rich material onto the surface of BD+48 740, since the RV measurements 
do not show any sign of it. 

Lithium could have been produced in the interior o f BD+48 740, brought up to the stellar 



surface and survived, if a Cameron-Fowler mechanism ( Cameron & Fowled 1 1971 ) had been in 



place in this low mass RGB star. For that, high mixing rates between the surface and the Li-forming 
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regions are required. It h as been shown by seve ral works now that Li-rich stars can be found all 
along the RGB (see e.g. iLebzelter et al.l 1201 2|) thus excluding the Charbonnel & Balachandran 
hooch process, in which the extra mixing, and therefore the Li enrichment, should be associated 
with the RGB bump, when the molecular weight discontinuity created by the inward penetration 
of the convective envelope is erased by the advancement of the H-burning shell. Even allowing 
for uncertainties, the location of BD+48 740 on the HR diagram supports the conclusions of these 
works. 

Several other sources for extra m ixing have been suggested to explain Li-ric h giants, i.e. 

tidall y enforced enhanced extra mixing dDenissenkov & Herwig 2004) . thermohaline ( Charbonnel & Zahn 



20071) . and Magneto-Thermohaline Mixing (|Denissenkov et al 



20091) . In fact, they could be oper- 



ating all together in BD+48 740, if we assume a scenario in which a co mpanion has been tidally 
captured and engulf ed into t he ste llar s urface (see Villaver & Livioll2009 ). In particular, it has been 
recently shown by iGaraudl (|201l|) and iTheado & Vauclairl (|2012l) that the accretion of metal rich 
planetary material onto the surface of exoplanet host stars can increase the Li surface abundance 
of the target stars on timescales that depend on the mass and stellar structure and h ave been esti- 
mated to be short, of the order of a few Myr (see also iDenissenkov & Herwigll2004 ). Siess & Livio 
(I1999Q have shown that engulfment of substellar objects could also trigger Li enhancement and a 
mass-loss enhacement on the RGB. No trace of IR excess can be found in the 2MASS and WISE 
data for BD+48 740 (J-K vs. K-WISE[3.4, 4.6, 12 and 22 /an]). Archive IRAS observations, up- 
per limits mostly, s uggest however , that t he star is in the extended shell and relatively high lithium 
abundance phase of ISiess & Livid (|1999|) scenario. 



The projected rotation velocity of BD+48 740 is typical for its luminosity and effective tem- 
perature. The actual amount of stellar rotation velocity increase depends strongly on stellar struc- 
ture (envelope mass) and the (unknown) rate of planet evaporation. The observed rotation velocity 
of BD+48 740 neither contradicts nor supports the engulfment scenario. 



5. Conclusions 

We present evidence of a high lithium abundance, A(Li)=2.33 ± 0.04, in the RG star BD+48 
740. This value is well above the expected limit of A(Li)=1.5 for evolved stars. We also present 
multi-epoch, precise radial velocities for the star which, although sparse, show periodic variations 
which can be interpreted as a result of the Keplerian motion of a planetary mass companion with 
msinz = 1.6 My, on a highly eccentric orbit of e=0.67. We find no evidence of a stellar-mass 
companion to BD+48 740. 

Given the current evolutionary status of the star, its Li abundance and the planet's current 
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orbit, we discuss a possibility that BD+48 740 had a second planet in an innermost orbit that could 
have been engulfed by the star. This possibility, although not directly verifiable with the available 
data, allows a scenario in which both the high lithium abundance of the star and the planet's highly 
eccentric orbit could originate from a recent, violent dynamical event, likely representing a later 
stage in the evolution of the planetary system, which was similar to HAT P-13 while still on the 
MS. 

BD+48 740 is the first example of Li-overabundant giant star with a planetary mass compan- 
ion candidate. 
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